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Abstract 

Unstable ^'^C nuclei are produced as quasi-projectiles in ^^C+^^Mg collisions at 
E/A = 53 and 95 MeV. The decay of their short-lived states is studied with the 
INDRA multidetector array via multi-particle correlation functions. The obtained 
results show that heavy-ion collisions can be used as a tool to access spectroscopic 
information of unbound states in exotic nuclei, such as their energies and the relative 
importance of different sequential decay widths. 
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Energetic heavy-ion collisions have been extensively studied to extract infor- 
mation about the properties of nuclear matter under extreme conditions pLj. 
These studies have also shown that a large variety of isotopes is produced dur- 
ing the dynamical evolution of the reaction. Some of these isotopes, being far 
from the valley of stability, live temporarily and decay by particle emission. 
Their unbound states can then be isolated and studied by means of correlation 
techniques PISIHIISIIS] . As an example, proton- ''Be correlation functions have 
been recently measured in central Xe+Au collisions at E/A=50 MeV ^ to de- 
termine the spin of internal unbound states of the astrophysically important 
nucleus. In this respect, a collision between two heavy ions can be viewed not 
only as a tool to study nuclear dynamics but also as a laboratory to produce 
several nuclear species in one single experiment and study their spectroscopic 
properties. This aspect of heavy-ion collision experiments represents an im- 
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portant perspective to access information about the properties of very exotic 
nuclei. Among all exotic nuclear species that can be produced in nuclear reac- 
tions, can be considered as an especially interesting one [7II5] . An extended 
Fermionic molecular dynamics (FMD) approach [9] allows one to describe light 
nuclei with cluster and halo structures [TOj. Within the context of similar ap- 
proaches, Antisymmetrized molecular dynamics (AMD) calculations predict a 
molecular structures for the ground state of [11]. Exotic features, such as 
molecular states and clustering [T2|13lll4lll5fl6lll7lll8fl9f20j . have indeed at- 
tracted the interest of a large community and particle correlation analyses 
have been used to access such features experimentally |14lll5lll6|17lll8lll9|20j . 
However, previous studies of nuclei have hardly reconstructed those states 
lying above the particle emission threshold of 3.73 MeV [2T]. These states 
decay into a single final configuration constituted by two alpha particles and 
two protons {2a+2p). All intermediate states that can be formed starting from 
the decay of nuclei by charged particle emission (i.e. ^He, ^Li, ^Be, *Be, 
^B) are unstable. Reconstructing the decay paths of these states therefore 
requires the detection of all intermediate and final particles with high energy 
and angular resolution. Furthermore, a large solid angle detector coverage is 
required to completely detect all decay products, due to the excitation energy 
of these high lying states. 



In this article we study nuclei produced as excited quasi-projectiles (QP) 
in ^^C+^^Mg peripheral collisions at E/A=53 and 95 MeV. The experimen- 
tal data are collected with the INDRA multidetector array [22|23] . By means 
of three- and four-particle correlation functions we provide experimental evi- 
dence of sequential decay modes for excited states through the production 
of intermediate unstable ^B, ^Be and ^Be nuclei. Exploring the relative con- 
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tributions of these sequential decay processes to the total decay widths of 
^'^C states provides important spectroscopic information which is relevant to 
access branching ratios and spins. 

In the first section, the experimental setup and the event selection criteria are 
presented. The multi-particle correlation analysis techniques used to explore 
sequential decay modes are illustrated in section [2] by selecting the special case 
of ^^C states. The last section of this work is finally devoted to the study of 
and its unbound states decaying into 2a+2p . 



1 Experimental setup and event selection 

A 2 mg/cm^-thick ^^Mg target was bombarded by beams of ^^C at E/A=53 
and 95 MeV, with an intensity of a few 10'' pps, produced by the GANIL cy- 
clotrons. The data presented in this work were collected with the INDRA mul- 
tidetector array [22|[23] . The data acquisition trigger required an event multi- 
plicity greater than 2. This trigger condition strongly favors peripheral events. 
Fig. [T] shows the velocity distribution along the beam direction, V //labf^beam 
(normalized to the velocity of the beam, Vteam) vs the charge of the detected 
particles at E/A=53 (left panel) and 95 MeV (right panel). These distributions 
show that there are two sources of emission of charges Z>1 having parallel 
velocities centered around values of V //lab ~0.2 Vbeam and V //lab ~0.9 Vbeam, 
respectively. In order to select excited Carbon (Z=6) QP breakup events, we 
take particles with V//iab '>Nbeam.l'^ and we also require that the sum of their 
charges is equal to six units. Among all the selected channels of Carbon QP 
breakup, we focus on those corresponding to the 3-alpha decay of ^^C* (i.e. 
i^C* ^3a) and to the 2a-2p decay of ^"C* (i.e. ^°C* -^2a+2p). N niab I-beam 
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distributions for Z = 1-2 particles in the selected ^^C* and ^"C* breakup events 
are shown on Fig. O The discontinuity at \ / /iab=^beam/'^ is due to the selec- 
tion criteria. More details about isolating the QP decay products can be found 
on Refs. pll25] . 

Fig. E] shows distribution of E^, defined as the total kinetic energy in the center 
of mass of the 2a-2p (left panels) and So; (right panels) exit channels. These 
distributions display several peaks, more pronounced in the case of the 3a 
exit channel at lower incident energies. The position of these peaks does not 
depend on the incident energy. 



2 Multi-particle correlations 

In order to illustrate our analysis techniques, we first study the most peripheral 
events where excited ^^C QP are produced and are reconstructed by detecting 
three alpha particles with V//iab > Vbeam/2. We define the three-alpha particle 
correlation function, l-l-R(Efc), as 

l + RiE,) = J-'^f'\ . (1) 

^ uncorr\-'^k ) 

where the correlated yield spectrum, Ycorr{Ek), is constructed with 3 a parti- 
cles detected in the same event. This spectrum is sorted with respect to the 
quantity E^, already defined in Fig. [31 If a ^^C QP is produced at an excita- 
tion energy E*, it follows that Ek = E* + Q, where Q is the mass difference 
in the ^^C^ 3a channel. The uncorrelated yield spectrum, Yuncorr{Ek), can 
be constructed with different techniques. In previous works [5] a dedicated 
Monte-carlo simulation was used to determine the uncorrelated background 
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of multi-particle configurations. Here we rather use an event-mixing technique 
[26)127] , where the three uncorrelated alpha particles are taken from three dif- 
ferent events. The full dots on Fig.|3]show the resulting 3a correlation function 
in i2c+24Mg collisions at E/A=53 MeV. The Q- value for the 3a decay of ^^C 
is -7.27 MeV (see the decay scheme of ^^C in Refs. |28ll2T] ). The first peak 
observed on Fig. H] at Ek ~0.5 MeV corresponds to the 7.65 MeV state in ^^C. 
The second peak, centered at Ek ~2.2 MeV, can be associated to the 9.65 
MeV state, while the third bump at E^ —6 MeV results from the overlap of 
highly lying closely packed states. 

The widths of the observed resonances are much larger than the intrinsic 
widths of the correspondent ^^C states. For example, the FWHM (Full- Width- 
Half-Maximum) of the state at 9.65 MeV (E^ ~2.2 MeV) is observed to be 
about 1.0-1.1 MeV, to be compared to an intrinsic width of 34 keV for this 
state. This discrepancy is mainly due to the finite angular resolution of the 
apparatus. We have performed a Monte Carlo simulation in order to estimate 
the effect of the INDRA detector on the resolution of the reconstructed excited 
states of ^^C decaying into three alpha particles. The excitation energy of ^^C 
was sampled continuously between and 15 MeV. The simulated events were 
filtered through the geometry and detector response of INDRA. Fig. [5] shows 
the distribution of the three-body kinetic energy, E^in — rec, as it is recorded 
by INDRA in the simulation, vs the input kinetic energy, Ekin- The average 
recorded kinetic energy, E^in — rec, is slightly shifted towards higher values as 
compared to Ekin- Also, the width in the Ekin — fee distribution increases with 
increasing Ekin- This is best shown on the bottom panel of Fig. [5] displaying the 
FWHM width of the Ekin — rec distribution as a function of Ekin- At Ekin='i-'i 
MeV, which corresponds to the 9.65 MeV state, the FWHM is about 1.0 MeV, 
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close enough to the measured one. From this study, it is clear that the position 
of the resonances is accurate but their widths are not. 

Stimulated by the analysis technique described in Ref. [5], we modified the 
definition of the uncorrelated three particle yields used in the denominator of 
Eq. [Hby using a partial event mixing (PEM) technique. In particular, the open 
symbols on Fig. H] correspond to the correlation function obtained when the 
denominator is constructed with two a particles taken from the same event 
and the third one from a different event. The resonant peaks are still observed 
in the correlation function, but with a reduced magnitude. This reduction 
can be partly attributed to the presence of 2-body correlations associated to 
sequential decays of ^^C. Indeed, the decay of these ^^C states can proceed 
via the emission of ^Be+a pairs, with the very loosely bound ^Be nucleus 
subsequently decaying into two a particles. These second step '^Be decays are 
included in the numerator of Eq. [T] but not in the denominator, if the latter 
is evaluated with the standard event mixing technique. In contrast, when the 
PEM technique is used, two a particles in the denominator are still taken from 
the same event. In this case, some a — a correlations due to secondary ^Be 
decays are kept in the denominator of Eq. [H Then, the sequential ^Be decay 
contributions contained in the numerator of Eq. [1] are partially cancelled out 
by calculating the ratio, Ycorr/yuncorr, and the magnitude of the correlation 
function peaks is reduced. Therefore, the observed attenuation of the peaks 
magnitude represents, by itself, an indication of the existence of ^^C sequential 
decay, ^"^C^^Be+a ^ a + a + a. 

In order to test the validity of the method, we construct the correlation func- 
tion for where no resonances are expected. For example, we select the 
four-body p-d-^He-a exit channel events. The quasi-projectile could decay 
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by first emitting ^He and ^Be. The ^Be nuclei is unbound but its states are not 
expected to decay into p-d-^He p9|30] . Therefore no peaks are expected in the 
correlation function. Fig. E] shows the correlation function for the three parti- 
cles constructed with the standard event mixing technique (full dots) and the 
PEM technique (open symbols). All combinations are free of correlations, the 
correlation functions are fiat and average to unity as expected. The enhance- 
ment of the correlation function observed at low values is not significant, 
mostly due to the very low statistics data collected in this low energy region. 
The observed slight tendency for l+R(Efc) to be less than unity at Ej. <2 MeV 
is also due to the suppression of the yields of detected events at low E^ val- 
ues, as it suggested by our Monte Carlo simulations. Difficulties in measuring 
R(Efc) at low Ejt are strongly due to the low angular resolution of INDRA. 

3 Multi-particle correlations and excited states 

We now turn to the study of ^'^C states. Fig. [7] shows the four-particle 2a- 
2p correlation function, 1 + R{Ek), constructed by using both the standard 
(closed symbols) and partial event mixing techniques (open symbols) already 
described in the case of 3a correlation functions. The top (bottom) panel refers 
to an incident energy of E/A=53 MeV (E/A=95 MeV). The 2a-2p correlation 
function constructed with the PEM technique (full dots) shows two broad 
peaks. The Q-value of the 2p+2a decay of is -3.7 MeV (see the level 
scheme of in ref [2T|30] ). The first peak at Ek ~1.5 MeV can be associated 
to an overlap of states at E* ^5.2-5.4 and 6.6 MeV. The second peak at 
Ek ~5-5.5 MeV corresponds to unknown states around E* ^9 MeV ^Ij. The 
location of these peaks is independent of the beam energy. The PEM explores 
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the following decays sequences: A) ^'^C^^B+p with the nucleus further 
undergoing a decay ^B^ a+a+p (open crosses); B) ^'^C^^Be+a with the 
^Be nucleus further decaying into a+p+p (open squares); C) ^°C— s>^Be+p+p 
with the ^Be nucleus further decaying into a+a (open circles). 

All PEM correlation functions display an attenuation of the peak magnitude. 
In the case of the states of observed at ^5-5.5 MeV, the attenuation 
obtained with the PEM techniques is less pronounced as compared to the 
case of the lower lying state at Ejt ?5il.5 MeV. This observation suggests that 
direct four body decays {2p+2a) without passing through any intermediate 
unbound state are more likely for the states around E* Ri9 MeV than for 
the state at E* ^^5.2-5.4 MeV. Among the studied sequential decay modes, 
sequence A induces a peak magnitude attenuation of the order of 30% with 
respect to its original value (see open crosses). Decay sequences B and C seem 
to provide comparable contributions to the decay of the studied state, 
within statistical uncertainties, reducing the peak magnitude to about 35- 
40% of its original value. A slight systematic preference for the B sequence 
as compared to the C sequence is observed. In general, for both the first and 
second group of peaks we observe a preference of to decay through the 
production of ^B, as compared to a decay proceeding through the ^Be system. 

In order to confirm these conclusions about sequential decay modes, we 
study three-particle coincidence spectra for events with 1 MeV< E^ < 3.5 
MeV, corresponding to the first peak on Fig. [71 On Fig. [8] (top panels) we 
show a — a — pi and a — a — p2 total kinetic energy spectra, N{Ek), where 
the total kinetic energy, Ek, is calculated in the three-body reference frame. 
In particular, we choose the slowest (pi) or the fastest (^2) proton in the 
reference frame. Similar spectra are constructed using the ai — p — p and 
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0(2 — P — P coincidences (bottom panels on Fig. [8]), where ai and a2 particles 
are chosen using the same criteria as pi and p2- These spectra of Fig. [8] are 
expected to contain information about possible intermediate states in and 
®Be nuclei produced during decay sequences A and B described above. 
If ^°C* did not decay through a definite sequence, the N (Ek) distributions 
would show no signatures of nuclear resonant decays. In contrast, we observe 
a peak at Ek ~0.2 MeV in the a — a—pi spectrum, corresponding to the decay 
of ^Bg s (decay sequence A), and a peak at E^. f»1.4 MeV in the ai — p — p 
spectrum, corresponding to the decay of ^Be^.^. (decay sequence B). In the 
spectra a — a — p2 and a2 — p — p (right panels on Fig. [8]) we observe broad 
peaks possibly associated to the decay of higher ^B and ^Be excited states. 
These observations seem to confirm our conclusions deduced from the PEM 
analysis shown on Fig. [71 The results obtained with our analysis techniques 
provide an experimental evidence of sequential decay modes of unbound states 
of ^''C exotic nuclei. Furthermore, our estimates of contributions from these 
different sequential modes have been recently confirmed by a dedicated exper- 
iment performed with a high resolution detector array and by using a different 
analysis techniques [7|. 

In conclusion, the existence of sequential decay modes of states proceed- 
ing through the production of intermediate loosely bound ^Be, ^Be and ^B 
nuclei has been explored. The analysis method, based on pure and modified 
event mixing techniques, has the advantage of not containing possible model- 
dependencies induced by the use of user-constructed simulations to build the 
uncorrelated background. Furthermore, the modified event mixing techniques 
allow us to semi-quantitatively estimate the relative contributions from differ- 
ent decay channels. A preference for a decay path producing unstable ^B nuclei 
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is found for the states at 5.2-5.4 MeV in ^'^C. The obtained results are also 
in good agreement with those extracted with high resolution experiments [7]. 
Our results confirm that heavy-ion collisions at intermediate energies can be 
used as a spectroscopic tool to access simultaneously important spectroscopic 
properties of several unstable nuclei all produced in one single collision. This 
perspective encourages an extension of similar studies to collisions induced 
by radioactive beams where more exotic nuclear species are expected to be 
produced. 
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Fig. 1. Scatter plots of charge (Z) vs velocity in the longitudinal direction, Vj/iab, 
normalized to the beam velocity, Vheam- The left panel refers to E/A = 53 McV 
beam energy; the right panel refers to E/A = 95 MeV beam energy. The dashed 
line indicates the cut used to select projectile break-up events (see text). 
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Fig. 2. Distributions of parallel velocities, V//;„b, normalized to the beam velocity, 
Vproj- for the ^^C QP decay (bottom panel) and for the ^'^C QP decay (top panel). 
The case of E/A = 53 MeV beam energy is shown on the left panels and the case of 
E/A = 95 MeV beam energy is shown on the right panels. The arrows show where 
the event selection cut was applied. 
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Fig. 3. Efe distributions for the ^'^C QP channel (left panels) and ^^C QP channel 
(right panels). Data corresponding to E/A = 53 MeV (95 MeV) beam energies are 
shown on the top (bottom) panels. 
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E, (MeV) 

Fig. 4. Three-a correlation functions. Pull dots: l+R(Efe) constructed with the stan- 
dard event mixing technique. Open symbols: l+R(Efe) constructed with the PEM 
techniques (see text). An expanded view of the correlation function at <5 MeV 
is shown in the inset. Particles pairs within parentheses are taken from the same 
event. 
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Ekin (MeV) 

Fig. 5. Monte Carlo simulation of the breakup of ^^C into three alpha particles. The 
top panel shows the distribution of kinetic energies, Ef^^n ~ f^c, as they are recorded 
by INDRA in the simulation, as a function of input Ekm values. The bottom panel 
displays the corresponding FWHM width of the Ekin—rec distributions as a function 

of Ekin- 
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Fig. 6. ^He-^H-^H correlation functions from the ^He-^He-^H-^H QP decay channel. 
Full dots: l+R(Efe) constructed with the standard event mixing technique. Open 
symbols: l+R(Ej.) constructed with the PEM technique. Particles within parenthe- 
ses are taken from the same event. 
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Fig. 7. 2a-2'p correlation functions at E/A=53 MeV (top panel) and 95 MeV (bot- 
tom panel). Solid symbols: standard event mixing technique. Open symbols: PEM 
techniques with two a's (circles), two a's and one proton (crosses) and one a and 
two protons (squares) taken from the same event (see text for details). Expanded 

views of the correlation functions for Eiq^5 MeV are shown in the inset. 
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Fig. 8. Top panels: a — a — -p kinetic energy spectra constructed with the slowest 
(pi, left panel) and the fastest (p2, right panel) proton. Bottom panel: a — p — "p 
kinetic energy spectra constructed with the slowest (ai, left panel) and the fastest 
(a2, right panel) a particle. The data correspond to an incident energy of E/A=53 
MeV. 
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